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We use self-consistent quantum transport theory to investigate the influence of electron-electron 
interactions on interlayer transport in semiconductor electron bilayers in the absence of an external 
magnetic field. We conclude that, even though spontaneous pseudospin order does not occur at zero 
field, interaction-enhanced quasiparticle tunneling amplitudes and pseudospin transfer torques do 
alter tunneling I-V characteristics, and can lead to time-dependent response to a dc bias voltage. 



I. INTRODUCTION 



In this paper we address interlayer transport in sep- 
arately contacted nanometer length scale semiconductor 
bilayers, with a view toward the identification of possible 
interaction-induced collective transport effects. When in- 
teraction effects are neglected a nanoscale conductor al- 
ways reaches a steady state^ in which current increases 
smoothly with bias voltage. Nanoscale transport the- 
ory has at its heart the evaluation of the density-matrix 
of non-interacting electrons in contact with two or more 
reservoirs whose chemical potentials differ. This problem 
is efficiently solved using Green's function techniques, for 
example by using the non-equilibrium Green's function 
(NEGF) methodi. Real electrons interact, of course, and 
the free-fermion degrees of freedom which appear in this 
type of theory should always be thought of as Fermi 
liquid theory^ quasiparticles. The effective single- 
particle Hamiltonian therefore depends on the micro- 
scopic configuration of the system. In practice the quasi- 
particle Hamiltonian is often^^^ calculated from a self- 
consistent mean- field theory like Kohn-Sham density- 
functional-theory (DFT). DFT, spin-density functional 
theory, current-density functional theory, Hartree theory, 
and Hartree-Fock theory all function as useful fermion 
self-consistent-field theories. Since a bias voltage changes 
the system density matrix, it inevitably changes the 
quasiparticle Hamiltonian. Determination of the steady 
state density matrix therefore requires a self-consistent 
calculation. 

Self-consistency is included routinely in NEGF 
simulations^^ of transport at the Hartree theory level in 
nanoscale semiconductor systems and in simulationsiSrii 
of molecular transport. The collective behavior captured 
by self-consistency can sometimes change the current- 
voltage relationship in a qualitative way, leading to 
an I-V curve that is not smooth or even to circum- 
stances in which there is no steady state response to a 
time-independent bias voltage. Some of the most use- 
ful and interesting examples of this type of effect oc- 
cur in ferromagnetic metal spintronics. The quasipar- 
ticle Hamiltonian is a ferromagnetic metal has a large 
spin-splitting term which lowers the energy of quasipar- 



ticles whose spins are aligned with the magnetization 
(majority-spin quasiparticles) relative to those quasipar- 
ticles whose spins are aligned opposite to the magne- 
tization (minority-spin quasiparticles). When current 
flows in a ferromagnetic metal the magnetization di- 
rection is altered. The resulting change in the quasi- 
particle Hamiltonian^^ is responsible for the rich va- 
riety of so-called spin-transfer torque collective trans- 
port effects which occur in magnetic metals and semi- 
conductors. These spin- transfer torque^iii are often 
understood macroscopically^^ as the reaction counter- 
part of the torques which act on the quasiparticle spins 
that carry current through a non-collinear ferromag- 
net. Spin-transfer torques^^ can lead to discontinuous 
I-V curves and to oscillatory^^ or chaotic response to 
a time-independent bias voltage. Similar phenomena^^ 
occur in semiconductor bilayers for certain ranges of ex- 
ternal magnetic field over which the ground state has 
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interlayer phase coherence or (equiva- 



lently) exciton condensation. Indeed when bilayer exci- 
ton condensate ordered states are viewed as pseudospin 
ferromagnets, the spectacular transport anomalies^^^ — 
they exhibit have much in common with those of ferro- 
magnetic metals. 

In this article we consider semiconductor bilayers at 
zero magnetic field, using a pseudospin language^^'^^ in 
which top layer electrons are said to have pseudospin 
up (I t)) and bottom layer electrons are said to have 
pseudospin down (| |)). Although inter-layer transport 
in semiconductor bilayers has been studied extensively 
in the strong field quantum Hall regime, work on the 
zero field limit has been relatively sparse and has focused 
on studies of interlayer drag,^^^^ counter flow, and on 
speculations about possible broken symmetry states. ^^'^^ 
We concur with the consensus view that pseudospin fer- 
romagnetism is not expected in conduction band two- 
dimensional electron systems^^'^^ except possibly^^ at 
extremely low carrier densities. There are nevertheless 
pseudospin-dependent interaction contributions to the 
quasiparticle Hamiltonian. When a current flows the 
pseudospin orientation of transport electrons is altered, 
just as in metal spintronics, and some of the same phe- 
nomena can occur. The resulting change in the quasi- 
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FIG. 1. Device cross-section in x, z direction. Top and bottom 
Alo.9Gao.1As barrier thickness is 60 nm, GaAs quantum well 
thickness is 15 nm, and Alo.9Gao.1As barrier thickness is 1 
nm. A channel is 1.2 \im long and 7.5 /xm wide. An assumed 
electron density in 2 dimensional electron gas (2DEG) is 2.0 x 



particle Hamiltonian is responsible for a current-induced 
pseudospin-torque which alters the state of non-transport 
electrons well away from the Fermi energy. Indeed al- 
though the spin-splitting field appears spontaneously in 
ferromagnetic metals, spintronics phenomena usually de- 
pend on an interplay between the spontaneous exchange 
field and effective magnetic fields due to magnetic-dipole 
interactions, spin-orbit coupling, and external magnetic 
fields. In semiconductor bilayers the pseudospin external 
field is due to single-particle interlayer tunneling ampli- 
tude and is typically of the same order^^ as the interac- 
tion contribution to the pseudospin-splitting field. 

The paper is organized in the following manner. We 
start in Section II with a qualitative discussion which 
explains the role played in transport by exchange corre- 
lation enhancement of the quasiparticle interlayer tunnel- 
ing amplitude. We also briefly outline the methods and 
the approximations that we have employed in the illus- 
trative numerical calculations which are the main sub- 
ject of this paper. In Section III, we describe the re- 
sults of self-consistent quantum transport calculation for 
an ideal disorder-free system with two GaAs quantum 
wells separated by a small AIGaAs barrier. For balanced 
double-layer electron systems, the inter-layer current has 
a resonant peak near zero-bias. Because the current is 
a non-monotonic function of inter-layer bias voltage, it 
is not immediately obvious that it is possible to drive a 
current that is large enough to induce exchange related 
instabilities. The calculations in Section III demonstrate 
that interactions alter the conductance peak near zero 
bias which occurs in weakly disordered systems in the 
absence of a strong magnetic field, and that exchange 
interactions can destabilize transport steady states. Fi- 
nally, in Section IV, we summarize our results. 



II. BILAYER PSEUDOSPIN TORQUES 
A. Interacting Bilayer Model 

We begin this section by briefly describing the ap- 
proximations that we use to model semiconductor bi- 
layer tunneling I-V characteristics. We consider an 
Alo.9Gao.1As/GaAs bilayer heterostructure with 60 nm 
top and bottom Alo.9Gao.1As barriers which act to iso- 
lated the coupled quantum wells from electrostatic gates, 
as illustrated in Fig. [H The bilayer consists of two 
15nm deep GaAs quantum wells with an assumed 2DEG 
electron density of 2.0 x lO^^cm"^ separated by a Inm 
Alo.9Gao.1As barrier. The quantum wells are 1.2/im long 
and 7.5/im wide with the splitting between the sym- 
metric and antisymmetric states set to a small value 
^SAS = 2t = 2/ieV. We define the z-axis as the growth 
direction, the x-axis as the longitudinal (transport) direc- 
tion, and the y-axis as the direction across the transport 
channel as shown schematically in Fig. ([1]). We connect 
ideal contacts attached to the inputs and outputs of both 
layers. The contacts inject and extract current and enter 
into the Hamiltonian via appropriate self-energy terms. ^ 

We construct the system Hamiltonian from a model 
with a single band effective mass Hamiltonian for top 
and bottom layers and a phenomenological single-particle 
inter-layer tunneling term: 
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(1) 



The first term on the right hand side of Eq. ([T]) is 
the single-particle non-interacting term while the second 
term is a mean-field interaction term. In the second term 
in right hand side, represents the Pauli spin matri- 
ces in each of the three spatial directions [i = x^y^z, 
(g) represents the Kronecker product, and A is a pseu- 
dospin effective magnetic field which will be discussed 
in more detail later in this section. To explore inter- 
action physics in bilayer transport qualitatively, we use 
a local density approximation in which the interaction 
contribution to the quasiparticle Hamiltonian is propor- 
tional to the pseudospin-magnetization at each point in 
space. If we take the top layer as the pseudospin up 
state (It)) cind the bottom layer as the pseudospin down 
state (I I)), the single particle interlayer tunneling term 
contributes a pseudospin effective field H with magni- 
tude t = AsAs/"^ and direction x. In real spin fer- 
romagnetic systems, interactions between spin-polarized 
electrons lead to an effective magnetic field in the di- 
rection of spin-polarization. Bilayers with pseudospin- 
polarization due to tunneling have a similar interaction 
contribution to the quasiparticle Hamiltonian. Includ- 
ing both single-particle and many-body interaction con- 
tributions, the pseudospin effective field A term in the 
quasiparticle Hamiltonian is,^^'^^ 



A = (t + /7m^Jx + /7m^,^ 



(2) 
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where the pseudospin-magnetization nips is defined by 



mps = ^Tt[ppst]. 



(3) 



In Eq. (j3j) , T = cTa; , cr^ , (7;^ is the vector of Pauh spin ma- 
trices, and pps is the 2x2 Hermitian pseudospin density 
matrix which we define as, 
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(4) 



The diagonal terms of pseudospin density matrix 
(Ptt' Pu) electron densities of top and bottom lay- 

ers. In Eq. ([2j), we have dropped the exchange potential 
associated with the z component of pseudospin because 
it is dominated by the electric potential difference be- 
tween layers induced by the inter-layer bias voltage (See 
below). From the definition of Eq. ([1]), the pseudospin- 
magnetization of x^y^z direction is defined in terms of 
the density matrix as 



(5) 
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As a result, we may express the system Hamiltonian in 
terms of pseudospin field contributions. 
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Here A^ is the electric potential difference between the 
two layers which we evaluate in a Hartree approximation, 
disregarding its exchange contribution. The planar pseu- 
dospin angle which figures prominently in the discussion 
below is defined by 



= tan 
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(9) 



This angle corresponds physically to the phase difference 
between electrons in the two layers. 




FIG. 2. Schematic of the Fermi surfaces of the bilayer system 
containing populations of symmetric and antisymmetric state 
separated in energy by a single particle tunneling term, /S^sas- 



—X directions respectively. The majority and minority 
pseudospin states differ in energy at a given momentum 
by 2teff where 



teff = t + U- 



(10) 



The population difference between symmetric and anti- 
symmetric differences may be evaluated from the differ- 
ences in their Fermi radii illustrated in Fig. [2l 
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where is the density-of-states of a single layer. Com- 
bining Eq. (p!Q|) and Eq. (pT]) we can relate U to 5, the 
interaction enhancement factor for the interlayer tunnel- 
ing amplitude: 

= ^ . St. (12) 

The physics of S is similar to that responsible for the in- 
teraction enhancement of the Pauli susceptibility in met- 
als. According to microscopic theory a typical value 
for S is around 2. We choose to use 5* rather than U as 
a parameter in our calculations and therefore set 

C/=^^-:^. (13) 
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B. Enhanced Interlayer Tunneling 

With the quasiparticle Hamiltonian for our semicon- 
ductor bilayer defined, we now address the strength of the 
interlayer interactions present in the system. The inter- 
action parameter U in Eq. ([2]) is chosen so that the local 
density approximation for interlayer exchange reproduces 
a prescribed value for the interaction enhancement of the 
interlayer tunneling amplitude. In equilibrium the pseu- 
dospin magnetization will be oriented in the x direction 
and the quasiparticle will have either symmetric or anti- 
symmetric bilayer states with pseudospins in the x and 



C. Pseudospin Transfer Torques 

In the following section we report on simulations in 
which we drive an interlayer current by keeping the top 
left and top right contacts grounded and applying iden- 
tical interlayer voltages, V/att, at the bottom left and 
bottom right contacts. 

We choose this bias configuration so as to focus on in- 
terlayer currents that are relatively uniform. The trans- 
port properties depend only on the quasiparticle Hamil- 
tonian and on the chemical potentials in the leads. Be- 
cause the pseudospin effective field is the only term in 
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the quasiparticle Hamiltonian which does not conserve 
the z component of the pseudospin, it follows that every 
quasiparticle wavefunction in the system must satisfy^^ 
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where is the z component of the pseudospin current 
contribution from that orbital, i.e. the difference be- 
tween bottom and top layer number currents, and nips 
is the pseudospin magnetization of that orbital. For 
steady state transport, the quasiparticles satisfy time- 
independent Schroedinger equations so that, summing 
over all quasiparticle orbitals we find. 
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(15) 



In Eq. ([IS]), is the planar orientation of A. The first 
equality in Eq. (p!5|) follows from Eq.(j2j). The pseudospin 
orbitals do not align with the effective field they expe- 
rience because they must precess between layers as they 
transverse the sample. The realignment of transport or- 
bital pseudospin orientations alters the total pseudospin 
and therefore the interaction contribution to A. The 
change in u\^s x ^ due to transport currents is referred 
to here as the pseudospin transfer torque, in analogy with 
the terminology commonly found in metal spintronics. 
Integrating Eq. ([15]) across the sample from left to right 
and accounting for spin degeneracy, we find that 



(16) 



where A is the 2D layer area, the angle brackets denote 
a spatial average, and II and Ir are the currents fiow- 
ing from top to bottom at the left and right contacts. 
(The pseudospin current 'f fiows to the right on the right 
and is positive on the right side of the sample, but fiows 
to the left and is negative on the left side of the sam- 
ple.) If the bias voltage can drive an interlayer current 
/ larger than 4etA(m^g)/fi, it will no longer be possible 
to achieve a transport steady state. Under these circum- 
stances the interlayer current will oscillate in sign and 
the time-averaged current will be strongly reduced. In 
the next section we use a numerical simulation to assess 
the possibility of achieving currents of this size. 

A similar conclusion can be reached following a dif- 
ferent line of argument. The microscopic operator J 
describing net current fiowing top layer (t) to bottom 
layer (I) given by— 
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where k, <j are the momentum and spin indices respec- 
tively. Here we have introduced a common notatition 



^SAS = 2t for the pseudospin splitting between sym- 
metric and antisymmetric states in the absence of in- 
teractions and added a factor of 2 to account for spin- 
degeneracy. The pseudospin density operator m^^ in Eq. 
([T7|) is given by 



(18) 



As a result, current density fiowing in the given system 
may be written as 



— l^xyl sin^ 



(19) 



where \nixy\ = \J {rrixY + {^vY ^ is defined in Eq. 
([9]) and Jc is defined critical current density. Assuming 
the current fiow is uniform across the device, it is possible 
to obtain the critical current simply by multiplying Eq. 
([T9|) by the system area A to obtain 
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D. Simulation detail 
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Because we neglect disorder in this simulation we may 
write the Hamiltonian in the form of decoupled ID lon- 
gitudinal channels in the transport {x) direction, taking 
proper account of the eigenenergy of transverse {y) direc- 
tion motion^^. The calculation strategy follows a stan- 
dard self-consistent field procedure. Given the density 
matrix of the 2D system, we can evaluate the mean-field 
quasiparticle Hamiltonian. For the interlayer tunneling 
part of the Hamiltonian we use the local-density approxi- 
mation outlined in subsection A. The electrostatic poten- 
tial in each layer is calculated from the charge density in 
each of the layers by solving a 2D Poisson equation using 
an alternating direction implicit method^^ with appro- 
priate boundary conditions. The boundary conditions 
employed in this situation were hard wall boundaries on 
the top and sides of the simulation domain and Neumann 
boundaries at the points where current is injected to in- 
sure charge neutrality^. Given the mean-field quasipar- 
ticle Hamiltonian and voltages in the leads, we can solve 
for the steady state density matrix of the two-dimensional 
bilayer using the NEGF method with a real-space basis. 
The density-matrix obtained from the quantum trans- 
port calculation is updated at each state in the itera- 
tion process. The update density is then fed back into 
the Poisson solver and on-site potential is updated us- 
ing the Broyden method^^ to accelerate self-consistency. 
The effective interlayer tunneling amplitudes are also up- 
dated and the loop proceeds until a desired level of self- 
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FIG. 3. Plot of the height of the inter layer conductance as 
a function of the interlayer exchange enhancement S at in- 
ter layer bias Vt — Vb = 10 nV. In the inset, we plot of the 
height of the interlayer conductance normalized to the non- 
interacting conductance (S=l) as a function of S with the 
same x axis. We clearly see that the height of the interlayer 
conductance follows an dependence. The single particle 
tunneling amplitude t = A^as/S = 1 /j^eV. 
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FIG. 4. Plot of the differential conductance of the bilayer sys- 
tem at 5* = 4 (orange triangle) , 5* = 2 (blue circle) , and S = 1 
(opened rectangular) as a function of interlayer bias Vb — Vr 
at OK. In the inset, the same plot in negative interlayer bias 
with same y axis is plotted with x axis in log scale. In small 
bias window, almost constant interlayer transconductance is 
observed. After the interlayer bias reaches Vint ~ lO/iV, an 
abrupt drop in transconductance can be seen in the interact- 
ing electron cases of 5* = 2,4. We plot the current vs. bias 
voltage for each enhancement factor in the inset. 



consistency is achieved. The transport properties are cal- 
culated after self-consistency is achieved by applying the 
Landauer formula, i.e. by using 
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III. NUMERICAL RESULTS 

A. Linear response 

In Fig. (|3]), we plot the interlayer conductance at tem- 
perature T = as a function of the interlayer exchange 
enhancement, S by setting Vb = Vjnt and Vr = in 
Fig. (dl). The depend ence demonstrates that the 
conductance in our nanodevice is proportional to the 
square of the quasiparticle tunneling amplitude, as in 
bulk samples^^^^'^^'^^, and that the quasiparticle tun- 
neling amplitude is approximately uniformly enhanced 
even though the finite size system is not perfectly uni- 
form. The range of S used in this figure corresponds to 
the relatively modest enhancement factors that we expect 
in bilayer systems in which the individual layers have the 
same sign of mass. For systems in which the quasiparti- 
cle masses are opposite in the two layers we expect values 
of S that are significantly larger. Spontaneous interlayer 
coherence, which occurs in a magnetic field^^'^^^^ but is 
not expected in the absence of a field would be signaled 
by a divergence in S. 



The physics of the results illustrated in Fig. ([3]) can 
be understood qualitatively by ignoring the finite-size- 
related spatial inhomogeneities present in our simulations 
and considering the simpler case in which there is a single 
bottom-layer source contact and a single top-layer drain 
contact on opposite ends of the transport channels. The 
interlayer tunneling is then diagonal in transverse chan- 
nel, and the transmission probability from top layer to 
bottom layer in channel k is 



' interlayer 



sm 



(22) 



In Eq. (|22]) . Sk is the difference between the current di- 
rection wavevectors of the symmetric and antisymmetric 
states and L is the system length in the transport di- 
rection. We may simplify the expression in Eq. (|22]) 
by rewriting it in terms of the Fermi velocity, Vf in the 
transport-direction and making use of the small angle 
expansion of the sin function to obtain 
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where tg// is the quasiparticle tunneling amplitude pro- 
portional to the interlayer exchange enhancement S re- 
sulting in the power law dependence seen in Fig. (|3|). As 
this approximate expression suggests, we find that trans- 
port at low interlayer bias voltages is dominated by the 
highest energy transverse channel, which has the smallest 
transport direction Fermi velocity. 
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FIG. 5. Pseudospin polarization vs. current. The quantities 
are expressed in units of the value of rrix at zero bias (mo). 
The points in this plot correspond to the same points that 
appear in the inset of current vs. bias voltage in Fig.Q. For 
zero-current, the pseudospin is in the x direction and has a 
value close to mo = z^o^S't. rux decreases monotonically, ruz in- 
creases monotonically, and rriy increases nearly-monotonically 
before saturating at the largest field values. The arrow in the 
plot indicates the direction in which a bias voltage increases. 



B. Transport beyond linear response 

We have seen that at small bias voltages the inter- 
layer current is enhanced by interlayer exchange inter- 
actions. In Fig. dll , we compare our interlayer trans- 
port result with the 5* = 1 case. We see that all of the 
curves are sharply peaked near zero bias, as in bulk 2D 
to 2D tunneling^^^ — j^i^. In all cases the decrease and 
change in sign of the differential conductance at higher 
bias voltages is due to the build-up of a Hartree potential 
difference (A^) between the layers which moves the bi- 
layer away from its resonance condition. The peak near 
zero bias is sharper for the enhanced interlayer tunneling 
(filled triangle and circle in Fig. (j4])) cases. To get a 
clearer picture of > 1 transport properties in the non- 
linear regime, we examine the influence of bias voltage on 
steady-state pseudospin configurations. Initially, all or- 
bitals are aligned with the external pseudospin field, the 
inter-layer tunneling, and it follows from Eq. ([16]) that 
(j)ps = 0. The enhanced magnitude of this equilibrium 
pseudospin polarization is mo = v^St = voteff. When 
current flows between there must be a ^-component pseu- 
dospin, as we have explained previously, whose spatially 
averaged value is proportional to the current. If it is 
possible to drive a current that is larger than allowed by 
Eq.( [20]), it will no longer be possible to sustain a time- 
independent steady state. 

In Figs.© and (|6|) we plot the magnitudes of the x. 
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FIG. 6. Pseudospin polarization vs. bias voltage. The quan- 
tities are expressed in units of the value of rrix at zero bias 
(mo). The inset of this plot clearly shows that rUz increases 
monotonically with bias voltage which eventually causes the 
current to decrease. The \nixy \ decreases only slowly with bias 
voltage because of the compensating effects of pseudospin ro- 
tation toward the z direction and an increase in the overall 
pseudospin polarization. There is no steady state solution to 
the self-consistent transport equations beyond the largest bias 
voltage plotted here, because the pseudospin orientation 
has reached 7r/2. 



^, and z directed pseudospin fields, along with \mxy\ = 

^Jm'^ + m^, evaluated at the center of the device, as a 

function of current and bias voltage, respectively. The 
pseudospin densities are plotted in units of their equi- 
librium values. The fifteen data points for each of these 
curves correspond to the 15 data points in the current vs. 
bias potential plot in Fig. (|4]), so that the maximum cur- 
rent corresponds to a bias voltage of ~ 20/iV and the last 
data point corresponds to a bias voltage of ~ 100/iV. The 
first thing to notice in this plot is that the z pseudospin 
component increases monotonically with bias voltage, as 
a potential difference between layers builds up. This is 
the effect which eventually causes the current to begin to 
decrease. The x component of pseudospin increases very 
slowly with bias voltage in the regime where the current 
is increasing, but drops rapidly when current is decreas- 
ing. At the same time the y component of pseudospin 
evaluated at the device center rises steadily with current 
until the maximum current is reached and then remains 
approximately constant. When the two effects are com- 
bined \mxy\ first increases slowly with bias voltage and 
then decreases slightly more rapidly. The relatively weak 
dependence of Im^^^l on current can be understood as a 
competition between two effects. Because of the Coulomb 
potential which builds up and lifts the degeneracy be- 
tween states localized in opposite layers, the direction 
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FIG. 7. Quasiparticle layer-polarization parameter a as a 
function of interlayer bias. When the source-drain bias ex- 
ceeds around 40 iiV ^ the mismatch match between subband 
energy levels exceeds the width of the conductance peak, a 
increases rapidly, and the x — ^-plane pseudospin polarization 
and critical current decrease. The inset shows the same data 
in log scale to emphasize the low bias behavior. 

of the pseudospin tilts toward the z direction, decreas- 
ing the X — y plane component of each state. At the 
same time the total magnitude of the pseudospin field 
increases, increasing the pseudospin polarization. In a 
uniform system the two effects cancel when a i-direction 
pseudospin field is added to a uniform system. 

Our simulations suggest the following scenario for how 
the critical current might be reached in bilayers. The 
width of the linear response regime is limited by the life- 
time of Bloch states which is set by disorder in bulk sys- 
tems and in our finite-size system by the time for es- 
cape into the contacts. In the linear-response regime the 
current is enhanced by a factor of S'^ by inter-layer ex- 
change interactions. The maximum current which can 
be supported in the steady-state is however proportional 
to the bare inter-layer tunneling and to the x — y pseu- 
dospin polarization and is therefore enhanced only by 
a factor of S. From this comparison we can conclude 
that more strongly enhanced inter-layer tunneling quasi- 
particle tunneling amplitudes (larger S) increases the 
chances of reaching a critical current beyond which the 
system current response is dynamic. For the parame- 
ters of our simulation, the maximum current of around 
4nA is reached at a bias voltage of around 10/iV. At 
this bias voltage the average x — y plane angle of the 
pseudospin field is still less than 90°, and steady state 
response occurs. At the bias voltage increases further 
the total current decreases, but the pseudospins become 
strongly polarized in the z direction. We illustrate this 
behavior in Fig. ([7]) by plotting 

a = A,/VA2+A2 (24) 



VS. bias voltage, where is z directional pseudospin 
field and Ac is a constant beyond which charge imbal- 
ance of the system becomes significant^^. The magni- 
tude of the x — y direction pseudospin polarization conse- 
quently decreases. In our simulations the critical current 
decreases more rapidly than the current in this regime. 
As shown in Fig. ([5]), (j)ps approaches 7r/2 {rrix 0) at 
the largest bias voltages (~ 100/iV) for which we are able 
to obtain a steady-state solution of the non-equilibrium 
self-consistent equations. 



IV. DISCUSSIONS AND CONCLUSIONS 

For a balanced double quantum well system, the equi- 
librium electronic eigenstates have symmetric or antisym- 
metric bilayer wavefunctions. When the bilayer is de- 
scribed using a pseudospin language, the difference be- 
tween symmetric and antisymmetric populations corre- 
sponds to a x-direction pseudospin polarization When 
the bilayer system is connected to reservoirs that drive 
current between layers, it is easy to show that the non- 
equilibrium pseudospin polarization must tilt toward the 
^-direction. The total current that flows between layers 
is in fact simply related to the total ^-direction pseu- 
dospin polarization. These properties suggest the possi- 
bility that electron-electron interactions can qualitatively 
alter interlayer transport under some circumstances. In 
equilibrium interactions enhance the quasiparticle inter- 
layer tunneling amplitude, but not the total current that 
can be carried between layers for a given ^-direction pseu- 
dospin polarization. If the inter-layer quasiparticle cur- 
rent can be driven to a value that is larger than can 
be supported by the inter-layer tunneling amplitude, the 
self-consistent equations for the transport steady state 
have no solution and time-dependent current response is 
expected. This effect has not yet been observed, but is 
partially analogous to spin-transfer-torque oscillators in 
circuits containing magnetic metals. We have demon- 
strated by explicit calculation for a model bilayer system 
that it is in principle to induce this dynamic instability 
in semiconductor bilayers. 

The current-voltage relationship in semiconductor bi- 
layers is characterized by a sharp peak in dl/dV at small 
bias voltages, followed by a regime of negative differential 
conductance at larger bias voltages. In our simulations 
the pseudospin instability occurs in the regime of neg- 
ative differential conductance where dynamic responses 
might also occur simply due to normal electrical insta- 
bilities. It should be possible to distinguish these two 
effects experimentally, by varying the circuit resistance 
that is in series with the bilayer system. The interac- 
tion effect might also be more easily realized in bilayers 
in which the resonant interlayer tunneling conductance 
is broadened by intentionally adding disorder. 
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